arXiv: gr-qc/0112021v 1 11 Dec 2001 


PARAMETRIC GRAVITY WAVE DETECTOR 

G. GEMME 1 , A. CHINCARINI 1 , R. PARODI 1 , PH. BERNARD 2 and E. PICASSO 3 
l Istituto Nazionale di Fisica Nucleare, via Dodecaneso, 33, 1-161^6, Genova, Italy 
2 CERN, CH-1211, Geneva 23, Switzerland 
3 Scuola Normale Superiore, piazza dei Cavalieri, 1, 1-56126, Pisa, Italy 

Since 1978 superconducting coupled cavities have been proposed as a sensitive detector of 
gravitational waves. The interaction of the gravitational wave with the cavity walls, and the 
resulting motion, induces the transition of some energy from an initially excited cavity mode 
to an empty one. The energy transfer is maximum when the frequency of the wave is equal to 
the frequency difference of the two cavity modes. In 1984 Reece, Reiner and Melissinos built 
a detector of the type proposed, and used it as a transducer of harmonic mechanical motion, 
achieving a sensitivity to fractional deformations of the order Sx/x « 10~ 18 . In this paper the 
working principles of the detector are discussed and the last experimental results summarized. 
New ideas for the development of a realistic gravitational waves detector are considered; the 
outline of a possible detector design and its expected sensitivity are also shown. 


1 Introduction 

In this paper we shall discuss the mechanism of the interaction of a gravitational wave with 
a detector based on two coupled electromagnetic cavities. In previous works this issue was 
discussed using the concept of a dielectric tensor associated with the gravitational waveS The 
interaction was analyzed in the reference frame where the resonator walls were at rest even in 
presence of a gravitational perturbation. Here we shall analyze the interaction in the proper 
reference frame attached to the detector and we shall therefore consider both the coupling 
between the wave and the mechanical structure of the detector and the perturbation induced 
on the field stored inside the resonator due to the time-varying boundary conditions. 

The proposed detector exploits the energy transfer induced by the gravitational wave between 
two levels of an electromagnetic resonator, whose frequencies and lu 2 are both much larger 
than the characteristic angular frequency II of the g.w. and satisfy the resonance condition 
C 02 — uq = fi. The interaction between the g.w. and the detector is characterized by a transfer 
of energy and of angular momentum. Since the elicity of the g.w. (i.e. the angular momentum 
along the direction of propagation) is 2, it can induce a transition between the two levels provided 
their angular momenta differ by 2; this can be achieved by putting the two cavities at right angle 
or by a suitable polarization of the electromagnetic field axis inside the resonator. In the scheme 
suggested by Bernard et al. the two levels are obtained by coupling two identical high frequency 
cavities® The angular frequency lo\ is the frequency of the level symmetrical in the fields of 
the two cavities, and 0 J 2 is that of the antisymmetrical one. The frequency difference between 
the symmetric and the antisymmetric level is determined by the coupling and can be adjusted 
by a careful resonator design. Since the detector sensitivity is proportional to the square of the 
resonator quality factor, superconducting cavities must be used for maximum sensitivity. 



The power transfer between the levels of a resonator made up of two pill-box cavities, 
mounted end-to-end and coupled by a small circular aperture in their common end wall, was 
checked in a series of experiments by Melissinos et ah, where the perturbation of the resonator 
volume was induced by a piezoelectric crystal!’! 

Recently the experiment was repeated by our group with an improved experimental set-upl 
we obtained an order of magnitude sensitivity to fractional deformations of the resonator length 
as small as 51/1 ~ 10 -20 Hz -1 / 2 . 

In this paper we shall discuss the detector’s working principles and briefly review the last 
experimental results obtained by our group on the first prototype. Finally a possible detector 
design, based on two coupled spherical cavities, is discussed and its expected sensitivity is shown. 


2 Fundamental principles 

When the resonator’s boundary is deformed by an external force the local displacement vec¬ 
tor, u(r,t), can be expressed as a superposition of the mechanical undamped normal modes 
{£ Q (r)}: u(r,t ) = Yl a Qa(t)£a{^)- Qa{t) is the generalized coordinate of the a t h mode, obeying 
the dynamical equation of motion: 


Qa(t) + g a (t) + ulq a (t) = 

^ la. 

and the modes are normalized according to the relation 


( 1 ) 


'Vol 


p(x)£ a (x) ■ ip(x) dV = M a 5 a p 


( 2 ) 


being p{x) the mass density and M a the reduced mass of the a th mode. For a homogeneous 
system M a = M, where M is the mass of the system. 

In eq. [l] an empirical damping term, proportional to the velocity, has been included. 
fait) is the generalized force, given by 


fait ) = / fir, t) ■ iair) dV 

JVol 


( 3 ) 


where fir, t) is the external force density acting on the system. 

For a plane g.w travelling along the z axis the force density, in the proper reference frame 
attached to the detector, has the form: 


fix,t) = p{x) (h\ 


\x + h\y, h\x — h\y , 0 


( 4 ) 


where is the adimensional amplitude of the wave, and h\ = —h\, h\ = h\. 

To study the mechanism of the energy transfer between the two levels of an electromagnetic 
resonator perturbed by a gravitational wave we shall make use of the fact that the electromag¬ 
netic field inside the resonator can be expanded over the fields of the normalized, orthogonal 
normal modes {E n if)} and {H n ir )}0 
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For simplicity we shall assume that in the frequency range of interest only two e.m. modes 
give a significant contribution (n = 1,2), and that the external force couples strongly only to 
one mechanical mode (a = m). If we now introduce a perturbation of the resonator boundary 
and assume that the perturbation is small i.e. that we can still expand the fields inside the 



perturbed volume over the normal modes of the unperturbed resonator, we obtain the following 
set of equations for the magnetic field expansion coefficients: 


Hi(t) + w 2 (i) + wiWi(t) = -uiq m (t) (CE«i(t) + C&H 2 (t)) (6) 

yi 

W 2 (*) + %-Mt) + w|W 2 (t) = -U 2 q m (t) (C%Hi(t) + C^H 2 {t)) (7) 

y 2 

) ) Wm9m(} ~ "IT + “DF - (8) 

where we have defined the time-dependent expansion coefficients as £i(t) = y/e^ f v E ■ E^dV , 
'Hi(t) = y/JIo f v H ■ Hi dV and the coupling coefficients as: 


C% = [ {ui&i • H, - UjEi • Ej) £ m • dS (9) 

J s 

The integral in eq. || is made on the unperturbed surface of the resonator. 

The electromagnetic quality factor Q n = G n /R s , takes into account the dissipation arising 
from the finite conductivity of the walls. R s is the material-dependent surface resistance of the 
walls, and the geometric factor G n of the nth e.m. mode is given by: 


G n — ui n p,Q 


JyH 2 dV 

J S H 2 n dS 


( 10 ) 


In the following we shall assume Q\ ~ Q 2 = Q. 

The term f^, in the r.h.s. of eq. |], describes the deformation of the walls induced by the 
stored e.m. fields, i.e. a back-action effect of the e.m. field on the detector’s boundary: it is 
well known that in a resonant cavity the stored magnetic field interacts with the rf wall current, 
resulting in a Lorenz force which causes a deformation of the cavity shapeQu The radiation 
pressure is given by: 

p L{r,t) = ^ (ji 0 H(r,t ) 2 - e 0 E(r,tf) (11) 

Expanding again the fields E and H in terms of the normal modes (eq. [|) we get: 


Pdr-t) = 

\ (HdtfHdrf + H 2 {t) 2 H 2 (r) 2 + 2H l {t)H 2 {t)H 1 {f) ■ H 2 (r)) ~ 

\ (Si(t) 2 E 1 (f) 2 + £ 2 (t) 2 E 2 (f) 2 + 2£ 1 (t)£ 2 (t)E 1 (r) ■ E 2 (f)) (12) 

Since H\{t) ~ exp(jcuit), H 2 (t) ~ exp(ju; 2 t), £\(t) ~ exp(juj\t) and £ 2 {t) ~ exp(ju; 2 i) 
only the cross-product terms will give a significant contribution at the resonance frequency 
= lo 2 — The other, rapidly oscillating terms, will just give an average deformation of the 
detector’s walls, determining a static frequency shift of the resonant modes, which can easily be 
compensated by an external tuning deviceU 

The generalized back-action force (cfr. eqs. ||| and |8j), acting on the mechanical mode 
of the structure, will be given by: 

rba _ 

Jm 

(H, ■ H„j £,. iS- if,(t)&(t) / (fT, • F 2 ) ( m -dS = 
l^W 2 (t)Wi(t) - lBS£ 2 (t)ft(t) (13) 





where we have introduced the coefficients A™[ = f s ^Hi ■ im'dS and £>21 = f s (Ei • E ' 2 j 
dS. 

The analysis of the system of differential equations ||-[8| can be simplified if we neglect the 
small perturbation on the initially excited e.m. mode (say mode 1), just taking into account 
the effects on the initially empty mode. Furthermore we shall consider the coupling between 
two TE modes of a resonator: for these modes we have vanishing electric field on the resonator 
surface and, as can readily be calculated, = 0 and A 21 ~ 

With this assumptions we can recast the coupled system of equations in the following form: 

E-2 (t) + —'H2(t) + w|?f 2 (t) = (14) 

Qm(t) + q m (t) + + I_^_H 2 (t)Hi(t) (15) 

The solution of quations [Tf] is straightforward when the back-action term is switched 
off in eq. [l5| and if we assume TLi(t) = 3fJ(Aexp(iwii)). In this case we have the following 
asymptotic solution (in the frequency domain): 

_ fM/M 

qm[ ’ (-u> 2 +u>l + iuu m /Q m ) [ ’ 

and 


H 2 M = 

1 V 2 ACR _ fmjw-uJ/M _ 

2 — u> 2 + lo\ + IUJUJ 2 /Q — (w - iOl) 2 + iO + iuimiu — OJl)/Q m 

1 U2AC2I _ fm{w + UJl)/M _ 

2 -iO 2 + io\ +IUJUJ 2 /Q -(iO + U>i) 2 + i 0 2 n + + L >1 )/Qm 


The field amplitude will be maximum when ui ~ ui 2 - If we design our detector so that 
u >2 — u}\ = io m we have: 


^ 2 (^ 2 ) 


1 trim fm(Wm) QQm 

2 ° 21 M u 2 u 2 m 


(18) 


Eq. |18| shows that the amplitude of field in the initially empty mode is proportional to the 
amplitude of the field in the excited mode A, and to the electromagnetic (Q) and mechanic 
(Qm) quality factor of the system. 

A detailed dicussion of the derivation of equations [6| |8] and of the noise sources affecting the 
performance of the detector will be given elsewhere. 


3 Experimental results 

The electromagnetic properties of a prototype detector, made up of two pill-box cavities, 
mounted end-to-end, and coupled trough an iris on the axis, were measured in a vertical cryo¬ 
stat after careful tuning of the two cells frequencies. In order to get maximum sensitivity we 
need to have two identical coupled e.m. resonators, or, in other words, a flat field distribution 
between the two cavities. The symmetric mode frequency was measured at 3.03 GHz and the 
mode separation was 1.38 MHz. 

In order to suppress the noise coming from the symmetric mode at the detection frequency, 
the transmission detection scheme, with two magic-tees, was used, as described elsewhercB 

In figure ||| the signal from the A port of the output magic-tee is shown for an input power 
Pi = 1 W and no adjustments made on the phase and amplitudes of the rf signal entering and 
leaving the cavity. The overall attenuation of the symmetric mode is R ~ —48 dB. 











Figure 1: Superconducting cavity mounted on the test cryostat. 
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Figure 2: Transmission of the symmetric mode (no optimization) measured at the A port of the output magic-tee. 


After balancing the arms of the two magic-tees in order to launch the symmetrical mode at 
the cavity input and to pick up the antisymmetrical one at the cavity output, with 1 W (30 dBm) 
of power at the X port of the first magic-tee, 6.3 x 10 -15 W (-112 dBm) were detected at the A 
port of the second one, giving an overall attenuation of the symmetric mode of R ~ —140 dB 
(see figure |3|). At a detection frequency of VL/2it ~ 1 MHz the sensitivity of the system is quite 
independent from the value of R, because of the high cavity Q. Nevertheless for lower frequencies, 
in a range H < 10 kHz, where astrophysical sources of gravitational waves are expected to exist, 
this noise source can become dominant and the achieved rejection is fundamental in order to 
pursue the design of a working g.w. detector in the 1-10 kHz frequency range. 

The cavity loaded quality factor was Ql = 1 x 10 9 at 1.8 K, and the energy stored in the 
cavity with 10 W input power was approximately 1.8 J (limited by the maximum power delivered 
by the rf amplifier), with both the input and output ports critically coupled (/?i « /3 2 ~ 1). 

To excite the antisymmetric mode a piezoelectric crystal (Physik Instrument PIC 140, with 
longitudinal piezoelectric coefficient /t( = 2 x 10 -1 ° m/V) was fixed to one cavity wall. The 
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Figure 3: Transmission of the symmetric mode in the optimized system, measured at the A port of the output 
magic-tee. Measurement taken with 1 kHz resolution bandwidth. 

driving signal to the crystal was provided by a synthesized oscillator with a power output in 
the range 2-20 mW (3-13 dBm). The oscillator output was further attenuated to reduce the 
voltage applied to the piezo by a series of fixed attenuators and a variable attenuator (10 dB 
step). The oscillator frequency was carefully tuned to maximize the energy transfer between the 
cavity modes. 

The signal emerging from the A port of the output magic-tee was amplified by the LNA (48 
dB gain) and fed into a spectrum analyzer. In figure |] an example of the parametric conversion 
process is shown. 

The minimum detected noise signal level at the antisymmetric mode frequency, with no 
excitation coming from the piezo, was P ou t(uj 2 ) = 5 x 10 -19 W in a bandwidth 5f = 100 Hz, 
giving a noise power spectral density V ou t( UJ 2 ) = 5 x 10 -21 W/Hz; the main contribution to this 
signal was the johnson noise of the rf amplifier used to amplify the signal picked from the A 
port of the output magic-tee. 

Taking into account the input and output coupling coefficients the sensitivity if the system 
is given by h rnin « 3 X 10~ 2 °(Hz)^ 1/2 

4 Future perspectives 

The second phase of the R&D program is focused on the development of a detector based on two 
spherical coupled cavities (see figure ||). In order to approach the interesting frequency range 
for g.w. detection, the mode splitting (i.e. the detection frequency) will be uj 2 — wi ~ 10 kHz. 
The internal radius of the spherical cavity will be r ~ 100 mm, corresponding to a frequency of 
the TE 011 mode uj ~ 2 GHz. The overall system mass and length will be M ~ 5 kg and L ~ 0.8 
m. The choice of these frequencies for the resonator and mode splitting will be also useful in 
order to test the feasibility of a detector working at ~ 200 MHz and at a detection frequency of 
« 1 KHz. 

A tuning cell, or a superconducting bellow, will be inserted in the coupling tube between the 
two cavities, allowing to tune the coupling strength (i.e. the detection frequency) in a narrow 
range around the design value. 

The choice of spherical cells depends on several factors: 

• From the point of view of the electromagnetic design the spherical cell has the highest 


















































Figure 4: Artistic view of the coupled spherical cavities with the central tuning cell. 


geometrical factor, and so the highest quality factor, for a given surface resistance. 

For the TE 011 mode of a sphere the geometric factor G has a value G ~ 850 17, while for a 
standard elliptical accelerating cavity the TMoio mode has a value of G ~ 25017. Looking at the 
best reported values of quality factor of accelerating cavities, which typically are in the range 
10 10 — 10 11 , we can extrapolate that the quality factor of the TEon mode of a spherical cavity 
can exceed Q ~ 10 11 . 

• From the mechanical point of view it is well know that a sphere has the highest interaction 
cross-section with a g.w. and that only a few mechanical modes of the sphere do interact 
with a gravitational perturbation (the quadrupolar ones)! 

The mechanical design is highly simplified if the spherical geometry is used since the defor¬ 
mation of the sphere is given by the superposition of just one or two normal modes of vibration 
and thus can be easily modeled. In fact the proposed detector acts essentially as a standard 
g.w. resonant bar detector: the gravitational perturbation interacts with the mechanical struc¬ 
ture of the resonator, deforming it. The e.m. field stored inside the resonator is affected by the 
time-varying boundary conditions and a small quantity of energy is transferred from the initially 
excited e.m. mode to the initially empty one, provided the g.w. frequency equals the frequency 
difference of the two modes. A possible design of the detector makes use of both the mechanical 
resonance of the resonator structure, and the e.m. resonance. This can be accomplished if the 
detector is designed in order to have the mechanical mode frequency equal to the e.m. modes 
frequency difference u> m = ui 2 — wi- In particular, for the detector designed to work in the 10 
kHz frequency range, the two lowest quadrupolar modes frequency will be approximately at 4 
and 17 kHz. The expected sensitivities of the detector for L 02 — uj\ = 4 kHz and u >2 — oj\ = 10 
kHz are shown in figures || and ||. In the calculation of the above curves the brownian motion 
contribution to the detector noise as well as the noise coming from the detection electronics has 
been taken into account. Note that, also when L 02 — ui\ / ui m (fig. ||j) the sensitivity of the 
system is fairly good. | 

“Actually the sensitivity of the system at 10 KHz is better than the sensitivity at 4 KHz. This is essentially 
due to the lower value of the brownian noise at higher frequency. 



h . = 1.2424e-023 



Figure 5: Calculated system sensitivity for a periodic source after 1 year integration time (w m — LO 2 Ct^l — 4 

kHz, Q = 10 11 , Qm = 5 x 10 3 , T =1.8 K). 

• The spherical cells can be esily deformed in order to remove the unwanted e.m. modes 
degeneracy and to induce the held polarization suitable for g.w. detection. 

The interaction between the stored e.m. held and the time-varying boundary conditions is 
not trivial and depends both on how the boundary is deformed by the external perturbation 
and on the spatial distribution of the helds inside the resonator. It has been calculated that 
the optimal held spatial distribution is with the held axis of the two cavities orthogonal to each 
other. Different spatial distributions (e.g. with the held axis along the resonators’ axis) give a 
smaller effect or no effect at all. 

• The spherical shape can be easily used to investigate whether the niobium-on- copper 
technique could be useful for the detector hnal design. 

The choice between bulk niobium or niobium-on-copper for the hnal detector design has not 
yet been made and is still under investigation. Both techniques present in principle advantages 
and drawbacks. A prototype of two coupled spherical cavities in bulk niobium will be built at 
CERN in 2002. A single cell, seamless, copper spherical cavity has been built at INFN-LNL by 
E. Palmieri and will be sputter coated at CERN. 

5 Conclusions 

A hrst prototype of the detector, made up of two pill-box cavities, mounted end-to- end, has 
been built and successfully tested. A detector based on two coupled spherical cavities is now 
being designed, and preliminar tests on nomal conducting prototypes are being made. The 
planned timeline is as follows: 

• In 2002 a bulk niobium detector (two spherical cavities, u = 2 GHz, Q = 10 kHz, hxed 
coupling) will be built at CERN; 

• In 2003 a variable coupling detector will be built and tested. 

If experimental results will be encouraging, by the end of 2003 a proposal for the realization 
of a g.w. detector, based on superconducting rf cavities will be made. 




h . = 2.0019e-024 



Figure 6: Calculated system sensitivity for a periodic source after 1 year integration time (u>m = 4 kHz, ui 2 —uii = 

10 kHz, Q = 10 11 , Qm = 5 x 10 3 , T = 1.8 K). 
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